The Bulk Viscosity of High- Temperature QCD 



Peter Arnold and Caglar Dogan 

Department of Physics, University of Virginia, 
Box 400714, Charlottesville, Virginia 22901, USA 

Guy D. Moore 

Department of Physics, McGill University, 
3600 University St., Montreal QC H3A 2T8, Canada 
(Dated: February 2, 2008) 

Abstract 

We compute the bulk viscosity Q of high-temperature QCD to leading order in powers of the 
running coupling a s (T). We find that it is negligible compared to shear viscosity rj for any a s that 
might reasonably be considered small. The physics of bulk viscosity in QCD is a bit different than 
in scalar (ft 4 theory. In particular, unlike in scalar theory, we find that an old, crude estimate of 
C — 15(| — Vs) 2 n gives the correct order of magnitude, where v s is the speed of sound. We also 
find that leading-log expansions of our result for £ are not accurate except at very small coupling. 



I. INTRODUCTION AND RESULTS 

Studies of collective flow at RH1C [1], particularly of elliptic flow, seem to be well de- 
scribed by nearly ideal hydrodynamics [2]. In fact, it has recently been claimed that these 
experiments prove that the quark-gluon plasma is the most nearly ideal fluid known, with 
a viscosity close to the conjectured lower bound on viscosities in any system [3, 4]. Such 
startling claims should be tested in any way we have available. This requires studying flow 
in heavy ion collisions using non-ideal hydrodynamics, that is, hydrodynamics including vis- 
cous effects [5] . It also would be valuable to know as much as possible about the theoretical 
expectations for viscosity in the quark-gluon plasma. 

In an ideal hydrodynamical treatment, the evolution of the plasma is determined by 
stress-energy conservation, d^T^ = 0, together with an equilibrium equation of state which 
relates the pressure to the energy density, P = P cq (e). This should work whenever the 
system is locally in equilibrium, which is the case in the limit of arbitrarily slowly varying 
flow velocity ui{x). When Ui{x) varies somewhat in space, the fluid will not be precisely 
in local equilibrium, which will modify the stress tensor. For slowly varying Ui(x), the 
corrections to the stress tensor Ty can be expanded in gradients of Ui. The leading order 
corrections are parametrized by two quantities, the shear viscosity r\ and the bulk viscosity 



where the expression is implicitly written in the instantaneous local rest frame (where T 0i = 



While we are really interested in the viscosities r] and ( of the quark-gluon plasma at 
temperatures T ~ 200 MeV, where the theory is far from weakly coupled, we only possess 
reliable tools for computing dynamical properties such as viscosities at weak coupling. 1 
Hopefully, extrapolating these results to strong coupling should give the right ballpark for 
the same quantities at moderately strong coupling, with uncertainties of perhaps a factor of 
a few. This motivates investigating 77 and ( at weak coupling. 

In a relativistic system, on dimensional grounds, both rj and ( must scale as rj, ( oc T 3 . A 
great deal of study has gone into the shear viscosity in QCD. It has been known for 20 years 
that the parametric behavior is rj ~ T 3 /(«g logfl/aj) [9, 10]; the leading coefficient was 
closely estimated in 1990 [11], and complete results now exist both at leading logarithmic 
order [12] and full leading order [13] in the QCD coupling a s . On the other hand, the 
calculation of the bulk viscosity has been completely neglected. To our knowledge, no paper 
in the literature even correctly states what power of a s it is proportional to. The purpose 
of this paper is to fill this gap, by computing the bulk viscosity in weakly coupled QCD at 
leading order in a s , using kinetic theory. We will only consider the case of vanishing (or 
negligible) chemical potential, [i = 0. 

In the next section, we will review the relevant physics of bulk viscosity, explaining why 



The lattice is a rigorous nonperturbative tool for studying thermodynamic properties of the quark-gluon 
plasma at strong coupling, but dynamical properties such as viscosities are hard to study on the lattice; 
see for instance, Refs. [6-8]. 
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Here Too refers to the heaviest zero-temperature (current) quark mass which is smaller than 
or of order the temperature T. We use the subscript zero to emphasize that to represents 
a zero-temperature mass and not a finite-temperature effective quasi-particle mass. We 
will see that the physics of bulk viscosity is much richer than that of shear viscosity. In 
particular, the conformal anomaly (i.e. scaling violations) and the corrections to quasi- 
particle dispersion relations due to interactions, both irrelevant for shear viscosity, are both 
essential pieces of physics for bulk viscosity. Particle number changing interactions also 
play a much larger role in bulk than in shear viscosity. These qualitative points have been 
anticipated by the pioneering work of Jeon and Yaffe [14, 15] on bulk viscosity in relativistic 
4 theory. However, we shall see later that there are some significant qualitative differences 
between bulk viscosity in 4 theory and in QCD. 

Section III will present the details of the calculation of bulk viscosity. Our discussion 
will at times be abbreviated, referring back to previous papers [12, 13], where much of the 
technology has already been presented. We will end with a discussion in section IV. However, 
for the impatient reader, we now present our main results. The coefficients, missing in 
Eq. (1.2), are presented in Fig. 1 and Fig. 2. Here, Nf is the number of flavors of quarks. In 
Fig. 1, all quark flavors are assumed to be massless (too <C & s T); in Fig. 2, all but one flavor 

1/2 

is assumed to be massless, with that one flavor's mass in the range a s T < too < a s '*T. A 
comparison of bulk viscosity and shear viscosity for three massless flavors is given in Fig. 3 
as a function of a s . The figure makes clear that neglecting bulk viscosity in favor of shear 
viscosity is actually quite a good approximation, not only at weak coupling but probably also 
at moderately strong, physically interesting couplings. Fig. 4 shows the ratio C/ a tv^ which 
at very small a s approaches a constant with corrections given by powers of (log(l/a! s )) -1 . 
The dashed line shows an old, crude estimate of the ratio of bulk to shear viscosity which 
will be discussed in Sec. IV. 

Throughout this paper, we will not attempt to project our leading-order results to cou- 
pling higher than a s ~ 0.3. In previous studies of diffusion constants [13], it was found that 
this is where different formulations of the effective kinetic theory, which were equivalent 
at leading-order in coupling, no longer agreed within a factor of 2, suggesting a complete 
breakdown of the perturbative treatment. 2 

II. PHYSICS OF BULK VISCOSITY 
A. Basic picture 

When a fluid is uniformly compressed, it leaves equilibrium. The energy density rises, 
but the pressure temporarily rises by more than what is predicted by the equation of state. 3 
Under uniform rarefaction, the pressure temporarily falls further than is predicted by the 



2 See in particular Fig. 4 of Ref. [13] at m^/T = 2.4 for 3-flavor QCD, which corresponds to a s = 0.3. 

3 That the pressure is higher during compression and lower during rarefaction is dictated by the second law 
of thermodynamics; if the pressure during compression were lower than in equilibrium, one could construct 
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FIG. 1: Bulk viscosity for massless QCD at several values of iVf, as a function of the coupling a s 
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Bulk viscosity due to nonzero quark mass 




FIG. 2: Bulk viscosity when it is dominated by a single quark flavor's mass, as a function of a s , 
for several values of N{. 



a perpetual motion machine of the second kind, which rapidly compressed a fluid (encountering a lower 
than thermodynamic pressure) and then slowly expanded a fluid (encountering the full thermodynamic 
pressure). This constraint, that ( is positive, is another way of seeing that £ must be proportional to 
the second power of the beta function (or other source of conformal invariance breaking), since the beta 
function can be of either sign. 
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FIG. 3: Shear versus bulk viscosity: rj/s and (/s (s the entropy density) as a function of a s , for 
iVf=3 QCD, neglecting quark masses. Bulk viscosity ( has been rescaled by a factor of 1000. 
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FIG. 4: The ratio C/ a tv f° r Nf=3 QCD, neglecting quark masses. The dashed line shows the crude 
estimate of (4.1) with (3.32). As a s — > (and leading-log approximations to the leading-order result 
become applicable), the ratio approaches the limit C/afi] — ► 0.973. 

fall in the density and the equation of state. The bulk viscosity quantifies the time integral 
of this extra shift in the pressure (per e- folding of expansion) . 

The change in pressure occurs because the fluid leaves equilibrium. The time scale for 
weakly coupled QCD to relax towards equilibrium is set by the rate Y ~ a^T log[l/a] for 
a typical particle (p ~ T) to randomize its momentum p. The faster the fluid equilibrates, 
the nearer to equilibrium it remains, so the smaller the shift in the pressure; therefore the 
viscosity should be proportional to e/Y ~ T 3 /a^ log[l/cv]. This naive estimate turns out to 
be parametrically correct for shear viscosity. 

However, it is wrong for bulk viscosity. The reason is that QCD (at high temperatures and 
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away from mass thresholds) is a nearly conformal theory and the bulk viscosity vanishes 
in a conformal theory for two reasons. First, uniform compression or rarefaction is the 
same as a dilatation transformation. In a conformal theory a dilatation transformation is a 
symmetry and so the fluid will not leave equilibrium. Therefore, ( must be proportional to 
the breaking of conformal invariance. 

Furthermore, in a conformal theory, even if the fluid is out of equilibrium, the pressure 
still does not deviate from the value given by the equation of state, which for a conformal 
theory is exactly P = e/3. This is just the tracelessness of the stress-energy tensor in a 
conformal theory. For instance, consider massless A0 4 theory, 

at finite A [with ( — h++) metric convention]. The Euler-Lagrange equation is, 

9 2 - ^ = , (2.2) 

and the stress tensor and its trace are 

T, v = d^d u <f> - Vlu ,C , T* = -(d<P) 2 - . (2.3) 

Multiplying the Euler-Lagrange equation by shows that Tjf vanishes up to a total deriva- 
tive, which averages to zero. This argument is only flawed because of the conformal anomaly, 
which arises because of the running of A with scale. The bulk viscosity coefficient will there- 
fore contain another power of the smallness of conformal invariance breaking 

Thus, in a nearly conformal theory, the bulk viscosity coefficient ( vanishes as the second 
power of the departure from conformality: one power because the departure from equilibrium 
is small, and another power because any departure from equilibrium has a small impact on 
the pressure. For massless QCD, conformal symmetry is broken by the running of the 
coupling, (3(a s ) ~ a 2 , and so the bulk viscosity is 

C ~ a s 2 log[l/a s ] X (a ' )2 ~ lo^fel ' (2 ' 4) 

as claimed before. The presence of quark masses also constitutes a breaking of conformal 
invariance provided m ^ T (otherwise there are no quarks in the thermal bath and the 
influence of the quark can be neglected). In this case the pressure deviates from the massless 
value by a relative amount ~ m^/T 2 , and 

T 3 /m|\ 2 m 4 

C ~ a s 2 log[l/« s ] X \ f^J ~ Ta s 2 log[l/a s ] ' (2 ' 5) 

For future reference, note that if one formally defines the pressure as P = Tu/3 and 
linearizes the hydrodynamic formula (1.1) about global equilibrium P = P cq (e), then the 
bulk viscosity parametrizes 

AP - v 2 s Ae = CV • u, (2.6) 

where v s is the velocity of sound, given by t> 2 = dP eq /de, and AP and Ae are the local 
deviations of pressure and energy density. 
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(a) (b) 
FIG. 5: Examples of (a) number conserving and (b) number changing processes in cj) 4 theory. 

B. Number changing processes: Comparison with cf> 4 theory 

There is one detail which this brief discussion has brushed over. Viscosities are typically 
determined by the slowest process which is required for relaxation to equilibrium. Certain 
departures from equilibrium can be very slow to equilibrate, due to the presence of almost- 
conserved quantities. 

For instance, when considering bulk viscosity in scalar A0 4 theory, Jeon found [15] that 
the total particle number equilibrates very slowly. The dominant process which randomizes 
momenta and determines the shear viscosity is shown in Fig. 5a, with rate T ~ A 2 T. In 
contrast, an example process which changes particle number, required for bulk viscosity, is 
shown in Fig. 5b. One might naively expect that the particle number changing rate from 
such processes is T ~ A 4 T, but this misses a soft enhancement. Number change primarily 
occurs between low energy excitations, where Bose stimulation increases the rate. The 
correct estimate is that the number of excitations relaxes at a rate T ~ A 3 T, but this is still 
parametrically small compared to the 2^2 scattering processes of Fig. 5a. This leads to 
the result ( oc AT 3 in 4 theory, up to logarithms [15]. In scalar theory, number- changing 
processes are the bottleneck for the relaxation to equilibrium characterized by bulk viscosity. 

The same does not occur in QCD (at vanishing chemical potential 4 ), however, because 
number changing processes are much more efficient in gauge theory. 5 The analog of Fig. 5 
is Fig. 6. Number change is relatively fast even among hard particles and occurs by 1^2 
splitting of a hard particle into two other hard particles during a small-angle collision, such as 
depicted by Fig. 6b. The small-angle collision rate is of order a s T, and the nearly collinear 
emission from such scatterings costs one extra factor of a s , giving a hard splitting rate 
T ~ OilT. For comparison, the rate for a hard particle to randomize its momentum through 
2^2 collisions is of order a 2 Tlog[l/a: s ], which is larger by a logarithm. One might then 
suppose that number change is still the bottleneck process for bulk viscosity (by a logarithm), 
that the relevant rate is therefore a 2 T rather than a 2 Tlog[l/a s ], and that therefore there 
should be no logarithm in the parametric formula Eq. (1.2) for (. This turns out not to 



4 For the Standard Model at finite baryon number chemical potential /i and finite quark mass m , the 
bulk viscosity would be very large. Compressing the system changes the temperature, which shifts how 
much of the baryon number is stored in each quark type, in equilibrium. The actual distribution of 
baryons between quark types approaches this equilibrium value only by weak interactions, leading to a 
bulk viscosity ( <~ ^ 2 mg/G|T 7 for [i < T in cases where this is the rate-limiting process. Note however 
that in the early universe yu/T ~ 1CT 9 is negligible, while in a heavy ion collision weak interactions can 
be neglected entirely and one should take the numbers for each quark type to be separately conserved. 

5 See, for example, Ref. [16] and the related discussion of photon Bremsstrahlung in Ref. [17]. 
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(a) (b) 
FIG. 6: Examples of (a) number conserving and (b) number changing processes in QCD. 

be the case, though, because 2^2 scattering processes exchange gluons between hard and 
soft momenta efficiently, and soft gluon number changing processes are efficient enough to 
prevent a particle number chemical potential from developing. In section HIE, we will 
show that, because of Bose stimulation enhancements for soft gluon emission from hard 
particles, the total rate of number- changing processes per particle is 0(a^ 2 T), which is 
parametrically faster than the log) rates discussed above. It is the latter, 0(a 2 T log) 

rates that will therefore be the bottleneck for equilibration and which will determine the 
QCD bulk viscosity. 

III. DETAILS OF THE CALCULATION 
A. Overview 

We now proceed with the details of the calculation of bulk viscosity. Our general approach 
and notation will follow [12]. To begin, note that, at weak coupling, there are long lived 
quasiparticles, and a kinetic theory treatment should be valid. The plasma is well described 
by a phase space density for each particle type, f(x,p), which can be expanded about a 
local equilibrium distribution f eq (x,p) as 

f(x,p,t) = f etl (x,p,t) + f 1 (x,p,t), 
f eq (x,p,t) = (exp^h^-p.u^))]^!)- 1 , (3.1) 

with P = T _1 and 7„ = (1 — u 2 )^ 1 / 2 . The departure from equilibrium is determined by the 
Boltzmann equation, 

^ + v p -V m f = -C\f\, (3.2) 

with C[f] the collision integral. Above, E p and 

v p = V p E p (3.3) 

are the energy and velocity of a particle with momentum p. 6 To study transport coefficients 
such as viscosity, we are interested in small departures from equilibrium in the hydrodynamic 



6 We use the general formula (3.3) [which can be understood as the group velocity of a wave packet] because 
we would like to make a general treatment of quasi-particles with some dispersion relation E p , and there 
is no need at this point to specialize, for example, to E 2 = p 2 +m 2 . 
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limit of slow variation in x and t. The left-hand side of (3.2) is explicitly small because of 
the derivatives, and so we may replace / by / eq there. The collision term must be expanded 
to first order, C[f] oc /i, noting that C[/ cq ] = by local detailed balance. 

It is convenient to analyze the problem in a local region, choosing an approximate rest 
frame where u(x) and the variation of /3(t) can be taken to be small. To first order in these 
small quantities, the derivatives appearing on the left-hand side of the Boltzmann equation 
are 

dMx,p,t) = -f (E p )[l±f (E p )} d,[P(t)(E p -p-u(x))} , (3.4) 

where fo is the Bose or Fermi distribution 

HE) = (e? E T I)" 1 - (3.5) 

The departure from equilibrium, in the case of bulk viscosity, arises because 

V-u = X^0. (3.6) 

In Sec. Ill B below, we will use the derivatives (3.4) and thermodynamic relations (in a 
treatment slightly generalizing that of Jeon and Yaffe [14]) to rewrite the left-hand side of 
the Boltzmann equation (3.2) in the form 



+ v p ■ V x f° = Afo(l±/o) X(x) q a (p) (3.7) 



dt v " eq 

for the case of isotropic compression or expansion, relevant to bulk viscosity. Here, a is 
a species label, and q a (p) represents how much a particle of type a and momentum p 
contributes to the AP — v^Ae of (2.6). 

The departure fi from local equilibrium, at linearized order, must also be proportional 
to X(x), and it is convenient to parametrize it as 

fi(x,p) = p 2 f (l±f )X(x) X a (\p\). (3.8) 

The function wm be a nontrivial function of the magnitude of momentum p = \p\, 

but (in the local rest frame) it is direction independent, because X is a scalar quantity. 
Defining 

S a (p) = -Tq a (p)f (l±f ) , (3.9) 
the Boltzmann equation can be written as 

S a (p) = [C X ]\p) , (3.10) 

with C the linearization of the collision integral, which we will give in Sec. HID. 

The bulk viscosity is then determined as the shift in the pressure induced by the departure 
from equilibrium x- As we shall discuss, this is an integral over p of x a times the same source 
S a already introduced, 

C = /3 3 E^/ -^S a (p)x a ( P )=(<S,x), (3.H) 

where v a is the multiplicity of species type a. The collision operator C is Hermitian under 
this inner product, and we may formally write, 

C = (S , C-'S) . (3.12) 

This can then be treated variationally, by the techniques presented in Refs. [12, 13]. 
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B. General formula for q a (p) 

It remains to determine q a , to establish the form of C, and to explain how the integral 
equations will be solved. We will treat q a first, since it is the most different from the problems 
already addressed in Refs. [12, 13]. The second term on the lefthand side of Eq. (3.2) is 

v P ■ VJ = /3/ (l±/ ) PiVpjVjUi (3.13) 

at linearized order. Specializing to isotropic compression or expansion, VjMj = (<5ij/3) V • u, 
this becomes 

v p ■ V x f = /3/o(l±/ ) V • u . (3.14) 

Unlike the case of shear viscosity, the term dtf is also nonzero; the compression or 
rarefaction of the fluid causes its density, and therefore its temperature, to change with 
time. By the chain rule, 

df dpdf dp djPEp) 

~dt ~ Tt^ ~ ~ It /o(1±/o) ^~ • (3 ' 15) 

Now, stress-energy conservation implies 

8^ = -»• d t e = -(e + P)V-u. (3.16) 
By standard thermodynamic relations, 

e + P = T^ (3.17) 

(recall that P = —F for a theory without chemical potentials), and by the chain rule, 

de dP de ,„ „ . 

-r = = v~ 2 d t P . 3.18 

dt dt dP s K ' 

Combining (3.16) thru (3.18), 

d JL = _ V 2 V . uT ^- = v 2 s V -uf3^, (3.19) 
dt dT s 1 d(3 ' v ; 

and since the dependence of P on t is through its (3 dependence, it follows that 

f t =PvlV.u. (3.20) 

Therefore, combining (3.14), (3.15), and (3.20), the full lefthand side of the Boltzmann 
equation is 

(d t + v p ■ V x ) f = /3/ (l±/o) V • u - vf-^^j . (3.21) 
Comparing to the definition (3.7) of q a (p), we determine 

= ^ - (3.22) 



10 



A nice property of this formula is that one can easily verify that the source vanishes in 
a conformal theory. In a conformal theory, the only dimensionful scale would be T, and so, 
by dimensional analysis, E p must have the form E p = pF(p/T) for some function F. Using 
(3.3), (3.22), and the conformal result v% = |, one would then find q a (p) = 0. 

Before finding explicit expressions for E p and in QCD, let us briefly discuss the result 
for q a (p). The p ■ v p term represents the change v p ■ V^/eq in the quasiparticle distribution 
function due to free propagation. For all other transport coefficients we have computed 
[12, 13], this type of change was the appropriate "source" in the Boltzmann equation, and 
the collision integral was to be equated with it. But here the "source" has nonvanishing 
energy, and energy is conserved. The collision integral has an exact zero mode associated 
with energy conservation; therefore collisions will not erase the change in /, but will re- 
distribute it until it looks like a shift in the temperature. The size of the temperature shift 
is fixed by energy conservation — that is, by the amount of energy the p ■ v p term introduces. 
Therefore, the true departure from equilibrium is the difference between this p ■ v p source 
term, and the temperature shift which carries the same total energy. This is the role of the 
second v^d((3E p )/d(3 term. In other words, considering the linearized collision operator C 
as an operator on the space of departures from equilibrium 5f, we must project the source 
p ■ v p into the subspace orthogonal to the zero mode of C (since the eigenvector of the zero 
mode is not actually a departure from equilibrium). 

As a check, we give a general demonstration in Appendix A that the source term deter- 
mined by (3.22) indeed carries no energy in the quasi-particle approximation we have used 
throughout. (One may also eschew generality and instead directly check with the explicit 
formulas for q a (p) given in the next section.) In the appendix, we also discuss in more detail 
why the quasi-particle approximation is justified for a leading-order calculation of the bulk 
viscosity. 

With this in mind, we can see why it is this same q a {p) which is relevant in determining 
the pressure shift due to the departure from equilibrium oc x a (p)- Naively, the extra pressure 
due to a departure from equilibrium fi(p) should be | f p • v p . However, a general 

shift in the equilibrium distribution function by j\ leads to a shift in the energy. Bulk 
viscosity involves the difference between the actual pressure, and the pressure determined 
by e and thermodynamics, P(e). Therefore, we must subtract off (dP/de)5e = v^5e, the shift 
in the pressure due to the extra energy density contributed by f\. That is precisely what 
the second term in Eq. (3.22) does. 7 

C. Specific formula for q a (p) 

Now we will determine in detail the form q a {p) takes in QCD at weak coupling a s C 1. 
For simplicity we will also take quark masses m <C T, though nothing in principle stops us 
from considering the case of quark masses mo ~ T. We will assume that there is at most 
one quark species with non- negligible quark mass, which we denote M . In this case, the 



7 This subtraction is technically unnecessary if one has already projected the source to be orthogonal to 
the zero mode, since then no shift in the energy would be produced. However, it is convenient, because 
it allows a symmetric treatment of the source and the pressure shift, as is manifested by the symmetric 
appearance of S in (3.12). 
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energy of a quasiparticle excitation of momentum p ^> gT, to first order in g 2 , is given by 



E 2 p = p 2 + rr& , 



m lo,a ^ark] = m 2 a + = m\ a + 



Cf^T 2 ___ 2 , ^ 2 T 2 

4 



rr&fehion] = (C A + JV f f P ) ^ = ^±_% 2 T 2 , (3.23) 

where m ,a is the mass of quark species a. The masses here are the corrections to the 
large p dispersion relations. 8 We have written these expressions in terms of group Casimirs 
so that they can be evaluated for a general group, and have also given the specialization 
to QCD, where the adjoint Casimir C A = 3, and the fermions are in a representation with 
Casimir C F = 4/3 and trace normalization t F = 1/2. Here N{ is the number of light Dirac 
fermions, or half the number of Weyl fermions. 
Using these expressions, q a becomes 

<i a = ^(y--sV + ^)) , (3-24) 

™l = mL,a " ^yaj > ( 3 - 25 ) 

which coincides with the results of Jeon and Yaffe [14] . 

The speed of sound can be determined by writing out the temperature dependence of the 
pressure. At order g 2 and Mq, the pressure of the QCD plasma is [20] 

p = (a + bg 2 [ii 2 =T 2 ])T 4 + cM 2 T 2 , 

2 2 

a = 180 {AdA + 7AW = ISO (32 + 2lNf) ' 

b = ^ (2d A C A + 5N { d F C F ) = ^ (48 + 2(W f ) , 

c = ^d F = ^ , (3.26) 

where d A = 8 and d F = 3 are the dimensions of the adjoint and fermion color representations. 
Using e = TdP/ dT — P, one finds 

9 dP dP/dT 1 26 „, 9x cM 2 , H . 

up to 0(g 5 ), 0(mlg 2 /T 2 ), and 0{rriQ/T 2 ) corrections. Here, 

2 _ fx 2 dg 2 g 4 f AN f t F -UC A \ g 4 f 2N,-33 \ 

^ g) = W 2 ]=mT 2 { 3 ) = lev {-^) (3 ' 28) 



Their relation to frequently-used zero- momentum masses are = = for gluons and 

m^, = 2m| for massless quarks, where inn is the Debye mass, mp is the analogous screening mass for 
quark exchange, and m p \ is the plasma frequency For further details, see, for example, Refs. [18, 19]. 
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is the beta function of QCD, which enters on taking the temperature dependence of g 2 into 
account 9 . Similarly, the quantities m 2 introduced earlier involve and (3(g 2 ), and are 



m 2 [quark] 
fn 2 [gluon] 



C V T 2 



m 



0,a 



(3(g 2 



(Ca + N t t F )r< 

6 



Collecting these results, and making the approximation 



1 ( p 

E~ r , 



1 



-v 2 (p 2 + rh 2 a ) ^((|- v 2 s )p 2 -\rh 2 a ) 



P 



valid for m 2 /T 2 < 1 and (3{g 2 ) < 1, we find 
q q ' a = \Av 2 \p + 



%Pi.9 2 )T 2 



rn, 



0,a 



P 



where 



Q 



\Av 2 



\Av 2 \p + 



Ca^a + N { cIfCf 



(3(g 2 ) T 2 



P 



-h(2d A C A + hd F C F N { )(3(g 2 ) + Q0d F M 2 /T 2 



(3.29) 



(3.30) 



(3.31a) 
(3.31b) 

(3.32) 



36tt 2 (4c/ a + 7d F Nf) 

Here, the Mq in (3.32) appears in the q for every species, but the mg a in the second term of 
(3.31a) only contributes to the possibly massive quark, for which mo, a = Mq. Note that, as 
promised, q a is proportional to the source of conformal invariance violation, either the beta 
function or the current quark mass. Because S oc q enters quadratically in Eq. (3.12), we see 
that C wm depend quadratically on the size of conformal invariance breaking, as claimed. 



D. Variational method and the collision integral 

It remains to specify the form of the collision integral, and to explain how it will be 
inverted to establish ( using Eq. (3.12). Since the details here are rather similar to the 
previous literature [12, 13], we will be somewhat brief in our discussion. First, define an 
inner product as in Eq. (3.11) (summation over species label and integration over momenta). 
Then the Boltzmann equation and bulk viscosity can be formulated variationally; define 

Q(x)=(x,s)-^{x,C X ), (3.33) 

and observe that SQ/5x = when \ satisfies the Boltzmann equation (3.10). Furthermore, 
the value (3.12) of ( is 2Q evaluated at this extremum: 

C = 2Q max . (3.34) 



We are implicitly taking Td/dT holding fx/T fixed. But we would get the same answer if we performed 
the derivative holding \i fixed; in writing g 2 [^ 2 — T 2 ] in Eq. (3.26), what we really mean is that there is 
explicit fi dependence in the g 4 term, of form B(3(g 2 )log(T 2 /fj, 2 )T 4 . Holding /x fixed, (3(g 2 ) arises from 
the T derivative of this logarithm. 
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A variational Ansatz for x will give a lower bound on the value of the extremum which will 
improve rapidly as the variational basis is increased. Therefore, we write a multi-parameter, 
linear Ansatz for xip)i i n terms of a set of basis functions. As we will discuss momentarily, 
xip) oc p at small momenta, and x grows no faster than ~ p 2 /T at large p . Therefore, we 
use a slight modification of the basis functions considered in [12] , 

0m(p) = ( T + p yr_ 2 » m = l...K. (3.35) 
The function x a {p) is then assumed to be of form, 

X a {p) = Y,^M- ( 3 - 36 ) 

m 

Within this variational Ansatz, the required inner products for Q are 



a,m 



X,C X ) = J2^ a lx b n , (3.37) 

where 

>(p)S a (p), 





abmn 


D ra 


= Va 

Jp 


mn 


= y a 

J P 



<j> m (p)[C ab <l>n}(p), (3.38) 

where C ab means the collision integral for species a when species b is out of equilibrium by 
the amount indicated by x h ■ Considering to be a rank NK column vector S and C^ n 
to be a NK x NK matrix, where N is the number of possibilities for the species index a, 
the bulk viscosity (3.12) is 

( = SC- 1 S. (3.39) 

In practice, N = 2 (quarks vs. gluons) if all quarks are massless, and N = 3 (massive quark 
vs. massless quarks vs. gluons) if one quark species is massive. 

The detailed form of C is given in Ref. [13], 10 which we summarize here for completeness: 

Ctn = ^£ / \Mt f {p,k-p\k')\\^f5^{P+K-P'-K') 

8 cdef J P k P' k ' 

xmp)f$(m±wm±f£(k')} 

x [0 m (p)5 ac +0 m (A;)^-0 m (p')5 ae -0 m (A;')5 a/ ] 
x [0„(p) ( 5 f ' c +0 n (A;)5 M -0 n (p / )5 6e -0„(A;')^ / ] 

+tE 47F / ^ , ^^71(p , ;p^)5(p , -p-^)/o c (p')[i±/o(p)][i±/o(^)] 

1 cde J ° 



X 



[<pm{p'W c -<Pm{p)5 ad -<Pm{k)5 ae ] [<t> n {p')5 bc -<t> n {p)5 M -4> n {k)5 be ] . (3.40) 
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See specifically Eqs. (2.22) and (2.23) of Ref. [13] with xt.j replaced by X a (p) to specialize to the isotropic 
(I = 0) angular dependence relevant to bulk viscosity, and then define C^ n as in (3.38) of this paper. 
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All factors of the number of degrees of freedom of each species are implicitly included in these 
sums. 11 The detailed expressions for the 2^->2 amplitude M. and the effective l<->2 splitting 
function 7 fill two appendices of Ref. [13] and will not be reproduced here. In treating the 
kinematics of these processes, we have neglected the masses of all external states, which is 
consistent with our approximation, m <C T. In principle there is no obstacle to treating 
the case m ~ T, but we have not done so, primarily out of laziness. 

Besides the difference in the source, which we have already stressed, the other difference 
between bulk and shear viscosity calculations is in the angular dependence of x i* 1 the 
collision integral. For shear viscosity, it was not <fi m (p) which appeared above, but 4> m {p) PiPj- 
(See, for example, Ref. [13] for a discussion in the conventions of this paper.) When suitably 
averaged over the indices ij, this led to angular factors of P2{cos 6 p k) in the cross-term 
between Xmip) an d Xn(k), for instance, where P2(x) = (3x 2 — 1)/2 is the second Legendre 
polynomial. Since bulk viscosity arises due to X = V • u, a scalar quantity, this angular 
dependence is absent. This makes the calculation of the collision integral somewhat simpler, 
but it does add two complications involving zero modes of the collision operator, to which 
we now turn. 

E. Zero modes of C 

The first term in the collision integral (3.40), corresponding to 2 <-> 2 processes, has 
two exact zero modes, corresponding to all x a (p) 00 1 an d all x a (p) °t Pi corresponding 
to particle number conservation and energy conservation, respectively. The second term, 
corresponding to collinear 1 <-> 2 processes, breaks particle number but still has the zero 
mode corresponding to energy conservation. Therefore, the collision matrix C will have a 
zero mode, and can potentially have a second approximate zero mode to the extent that 
the 2 <-> 2 term is larger than the 1 <-> 2 term. Since the collision integral must be inverted 
in evaluating Eq. (3.39), we must address the exact zero mode. We will see that making 
a leading-log expansion of bulk viscosity (if such is desired) requires treating the 2^2 
term as larger, by a logarithm, than the 1 « 2 term for p ~ T. In order to understand 
why number-changing processes are not a bottleneck for equilibration, and to understand 
expansions in [log(l/a s )] _1 , we will need to address the approximate zero mode as well. 
Both of these zero modes are specific to the case of isotropic xip)i an d neither is relevant 
to the analysis of other standard transport coefficients such as shear viscosity and flavor 
diffusion constants. 

The presence of an exact zero mode in the collision integral is not problematic, precisely 
because the source S carries precisely zero energy, and so is orthogonal to the zero eigen- 
vector. Therefore, our previous expressions should be understood as valid in the subspace 
orthogonal to the zero mode of C. In practice our basis of functions § a m are not restricted 
to this orthogonal subspace. But the collision integral can be rendered invertible without 
changing its behavior in the orthogonal subspace by adding a constant times the projection 



In the convention of Ref. [13], the sums (no averages) over all initial and final colors are included in |.M| 2 
and 7, and each of the indices cdef in the explicit sums above denote gluons vs. different flavors of quarks 
vs. different flavors of anti-quarks. 
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operator for the pure temperature fluctuation (the zero mode); 

C±^C a ± + \(v a J ^i? P m (p)/o a (l±/o a )) / ^^0n(A:)/o 6 (l±/o)) , 

for any positive A. This renders C invertible; and while C 1 is A dependent, C 1 S is not, 
since S has zero projection onto the modified direction. In our numerical evaluations we 
have checked explicitly that the determined value of ( has no sensitivity to the added value 
of A. 

Next, consider the possible approximate zero mode, x a (p) a constant, corresponding to a 
chemical potential for particle number. 12 First note that the constant value must be the same 
for fermionic and bosonic species, because the set of 2 <-> 2 processes includes fermionic pair 
annihilation to gluons, which contributes at leading logarithmic order. However, no elastic 
2^2 scattering processes will drive a common chemical potential for both quark and gluon 
number to zero. For the case of bulk viscosity in A0 4 theory it was found that this played 
a major role in setting the bulk viscosity [14]. In that theory, x oc 1 is an approximate zero 
mode of the full collision operator: (x|C%) for x — 1 is parametrically small compared to 
typical hard collision rates. 13 

However, for the bulk viscosity of QCD, this would-be zero mode actually plays no role: 
the expectation (x\Cx) for x = 1 is parametrically large rather than small compared to 
typical hard scattering rates. The reason is that, while number changing collinear processes 
(the second term in Eq. (3.40)) are subdominant to 2 <-> 2 processes at generic momenta, 
they are very fast at producing and destroying soft gluons. To see this, let us estimate the 
total rate for a hard particle to produce a soft gluon of momentum k by Bremsstrahlung. 
Combine (i) the 0(g 2 T) rate for small-angle scattering, as in Fig. 6a, times (ii) a factor of g 2 
for absorbing or emitting the additional gluon in Fig. 6b, times (iii) an initial or final state 
factor of f(k) or 1 ± f(k) for that gluon, and (iv) a momentum integral dk/k (responsible 
for the logarithmically large rate of soft bremsstrahlung emission in vacuum 14 ). f(k) ~ T/k 
for small k, and the result for the number changing rate T 1 ^ 1 1S then 15 

r tota. „ g , T J dk f{k) „ g , T 2 J dk (3 



By "particle number," we mean the sum of quark, anti-quark, and gluon numbers, not a difference like 
quark minus anti-quark number. 

This permits a simplification in 4 theory whereby one can avoid solving an integral equation and instead 
determine the leading-order bulk viscosity from a simple expectation value (x|Cx) for x cx 1 [14]. 
In vacuum, there is an additional logarithmic factor for bremsstrahlung from an ultra-relativistic particle — 
a collinear logarithm J d 2 k±/k 2 _ ~ \n(q/m), where q is the momentum-transfer in the underlying 2^2 
collision. In our case, the most frequent collisions are the small-angle ones, whose impact parameter is 
limited by Debye screening, and q ~ m ~ gT so that there is no collinear log enhancement. 
For k <C T, the Landau-Pomeranchuk-Migdal (LPM) effect plays no role in gluon emission, as discussed 
qualitatively in Sec. 5.2 of Ref. [19] . This is different than the case of soft photon emission due to the O(gT) 
thermal mass and scattering of the emitted gluon, either of which, for k <C T, causes loss of the multiple- 
collision coherence that produces the LPM effect. Here is a quick argument: For small k, the internal hard 
particle line in Fig. 6b is off-shell in energy by an amount of order AE = E p +k~E p — Ek ~ {m^+k\)/{2k). 
The formation time of the gluon is therefore of order (AE)^ 1 < k/m 2 <~ k/(g 2 T 2 ) 7 which is small compared 
to the time l/(g 2 T) between collisions when k <T. 
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The infrared divergence of the integral will be cut off by the effective thermal mass m ~ gT 
of the emitted gluon, so that 

f7 4 T 2 

ptotal _ y g 3 T _ (3^3) 

m 

As discussed earlier, this is parametrically faster than the 0(g 4 T log) rate to redistribute 
momenta between soft and hard particles, which is the bottleneck which determines bulk 
viscosity. The total rate r^^ 1 for creating or absorbing soft particles can therefore be taken 
as formally infinite for the purpose of a leading-order calculation of bulk viscosity. 

The same result can also be obtained, with some difficulty, from Eq. (3.40) of this paper 
together with eqs. (B1-B6) of Ref. [13], which determine the splitting functions 7. In par- 
ticular, the k <C p behavior of 7| g (p';j9fc) and 7q g (p';pA;) is 7 ~ (g A T)p 2 /k. If we substitute 
0m — <Pn — X = 1 i n to the 1 <-> 2 term in Eq. (3.40), the k integration for p ~ T then gives 
the linear divergence f dk/k 2 of (3.42). 

This means that a chemical potential is actually very rapidly thermalized by number 
changing processes. Any x(k) which falls more weakly than x{k) oc k at small k will lead to 
a divergent collision rate, meaning that such departures from equilibrium are so efficiently 
equilibrated that we need not consider them. Therefore we should restrict our Ansatz for x 
to only functions which are linear or higher powers of k in the soft region. This justifies our 
choice in Eq. (3.35). Within this subspace of functions Xi the 2^2 part of the collision 
integral has only one zero mode, that associated with energy conservation, which we have 
already discussed. Therefore the small a s behavior will indeed be ( oc a 2 T 3 / log[l/a s ], and 
one can perform an expansion in logarithms of the coupling if desired. 



F. Expansion in log[l/a s ] 

In Ref. [13] it was shown that an expansion in inverse powers of ln[l / a s ] works surprisingly 
well at small values of a s , if it is carried to next-to-leading order. As we have just seen, 
there is no obstacle to making a similar expansion here. We have done so, by following 
the procedure described in detail in Ref. [13], but we find that the expansion works much 
less well than in the case of shear viscosity and number diffusion. The reason is that the 
dominant physics in shear viscosity and number diffusion is angle change. The charge q a in 
that case is 1 or |p| times a nontrivial function of angle. The departure from equilibrium, 
x(p), has nontrivial angular dependence, but turns out to have very simple |p| dependence, 
so a one parameter Ansatz works very well. In a next-to-leading log treatment, one fixes the 
|p| dependence of x{p) using the leading-log part of the 2^2 processes and evaluates the 
collision integral using this fixed form of x{p)- This works because this functional form of 
x(p) is essentially correct, whatever collision processes are involved. 

For bulk viscosity, on the other hand, the charge q a changes sign as a function of the 
particle's momentum, as the 1/p and p terms in Eq. (3.31) change relative importance. 
The 1/p term is also larger for gluons than for quarks, due to their larger thermal masses; 
therefore, over most of the momentum range the quarks and gluons display opposite depar- 
tures from equilibrium. In QCD, the physics of bulk viscosity is primarily the physics of 
re-arranging the |p| dependence of particle distributions. This is what the number changing 
1 2 processes do best; so they play a much larger role in bulk viscosity than in shear. 
Indeed, unlike the case of shear, dropping the 2 <-> 2 processes and retaining only the number 
changing ones would still give a finite answer for ( — which in fact turns out to be within 
a factor of 2 of the leading-order answer over most of the range of a s we have considered. 
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FIG. 7: (Color online) Functional form of x(p) as a function of p, shown for quarks and gluons 
in massless Nf = 3 QCD at a s = .053 (or mv/T = 1). The three curves are the functional form 
using the leading-log 2^2 processes only, using the number changing processes only, and using 
all processes. 



QCD, N f = 


Leading-log A 


NLL ii* jT 







0.443 


7.14 


.151 


2 


0.638 


7.57 


.282 


3 


0.657 


7.77 


.286 


4 


0.650 


7.93 


.279 


5 


0.622 


8.06 


.263 


6 


0.577 


8.17 


.242 



TABLE I: Next-to-leading log bulk viscosity, ( = Aa?T 3 / ln[/U*/mrj], and £ calculated using only 
number changing collinear processes, Cl<-»2- All Nf quarks are taken to be massless. 



However, the detail of how they re-arrange the momentum distributions is different than for 
the elastic processes. Therefore the detailed p dependence of x(p) is quite different if only 
the leading-log 2^2 processes are considered, than if the full collision integral is used. We 
illustrate this difference in Fig. 7. This limits the range of validity of the expansion in logs 
to the regime where the 2 <-> 2 processes are much faster. But as we just said, the 1 <-> 2 
processes are more important to bulk viscosity than to shear, so this requires the logarithm 
actually to be large. Therefore the expansion in logs works poorly and should not be used 
in treating bulk viscosity. 

Another consequence of the quite nontrivial form of x(p) is that several basis functions 
must be used to get accurate numerical values of (. For instance, we need at least 5 basis 
functions to get 0.1% accuracy, something accomplished with two basis functions for shear 
viscosity. For this reason, the results presented in Fig. 1 and Fig. 2 are "only" good to about 
0.1%. 
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FIG. 8: (Color online) Bulk viscosity £> plotted against mo rather than a s for massless QCD. 
The dotted curves are the leading-log results; the dashed line on the right is the result neglecting 
everything but number changing collinear processes. 



For completeness, Table I lists the first two coefficients in an expansion in leading logs 
for massless QCD: 



Aa 2 T 3 

c- 



ml = 2mL[gluon] = (1 + lN ( )g 2 T 2 . 



(3.44) 



The table also contains the coefficient ( = Ca 2 T 3 we would obtain if we ignored all 2 <-> 2 
processes and considered only the number changing processes. To display the futility of 
using the next-to-leading log results, we compare them with the leading order results in 
Fig. 8. The failure of the next-to- leading log approximation by a factor of at least 1.5 for 
Nf = 3 or 6 at m D > 1.1 T corresponds to a s > 0.05. 



IV. DISCUSSION 



The physics of bulk viscosity in QCD is very interesting. The QCD plasma leaves equi- 
librium under compression or rarefaction only due to conformal symmetry breaking, and 
the bulk viscosity depends quadratically on the size of conformal symmetry violation (either 
through quark masses or the beta function). To find the departure from equilibrium one 
must include the forward-scattering corrections to dispersion relations, and must account 
carefully for the shift in the plasma temperature. The departure from equilibrium due to 
compression is of opposite sign for high and low momentum excitations, and of opposite sign 
at intermediate momenta p ~ ttT for quarks versus gluons. Collinear splitting processes ac- 
tually dominate the equilibration of the plasma except at very small coupling, although 
in the formal weak coupling limit, equilibration should be logarithmically dominated by 
2^2 scattering, annihilation, and Compton processes — with the proviso that soft gluon 
bremsstrahlung is also included, since it prevents the development of a chemical potential 
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for particle number. Putting this physics all together, one finds that the bulk viscosity 
at leading-log order (i.e. for exceptionally small coupling) is ( ~ a 2 T 3 / log[l/a: s ], with a 
leading coefficient of about 1. More practically, however, one can see from Fig. 1 that the 
complete result to leading order in powers of a s is roughly ( ~ 0.2 a 2 T 3 for any reasonable 
perturbative value of a s (0.02 < a s < 0.3). 

The practical import of bulk viscosity in QCD is very limited, however, in the regime 
where a perturbative treatment has any hope of applicability. We find that even for a s = 1/3, 
the bulk viscosity is hundreds of times smaller than the shear viscosity. In practice, this 
means that bulk viscosity can be neglected, whenever shear viscosity plays a role. For 
instance, the decay of a sound wave depends on the combination ( + 4i]/3; so one may drop 
the ( term to a very good approximation. The expanding QCD plasma in an ultra-high 
energy heavy ion collision is expected to be quite anisotropic, so shear viscosity again plays 
a role and bulk viscosity can be ignored. Similarly, while the expansion of the QCD plasma 
in the early universe should have been nearly isotropic, any flows in the presence of a phase 
interface — the only circumstances where nonequilibrium behavior may leave records in the 
early universe — are expected to be quite anisotropic, and again shear viscosity will be more 
important than bulk. 

Besides the quite elegant physics involved in the bulk viscosity of QCD, it also provides a 
nice example of the dangers of interpreting scalar field theory as a toy model for gauge theory, 
with A playing the role of g 2 . In massless A0 4 theory, Jeon and Yaffe showed [14, 15] that 
the shear viscosity behaves as i] ~ T 3 /A 2 , while bulk viscosity behaves as ( ~ AT 3 log 2 [l/A]. 
For shear viscosity, the scalar theory provides a successful toy model, missing only the 
logarithmic dependence. For bulk viscosity, although some of the physics is the same, scalar 
field theory is a misleading guide to gauge theory, getting even the power of the coupling 
wrong. The difference arises because number-changing processes in scalar theory are slow 
compared to processes which redistribute hard momenta (rate A 3 T vs. A 2 T); in QCD, they 
are fast (a^ 2 T vs. a 2 T). 

One consequence of slow particle number changing rates for scalar theory, observed by 
Jeon and Yaffe, was that the bulk viscosity did not match the crude estimate 

1577(1 -v 2 ) 2 (4.1) 

that had previously been made for scalar theory in the literature [23]. 16 (This same relation 
was found earlier by Weinberg for a photon gas coupled to hot matter [24].) However, these 
same estimates turn out to be parametrically correct for QCD, reproducing (1.2). In QCD, 
the bottleneck rate is the same for both shear and bulk viscosity, 

vl-\ = 0{P{g 2 ))+0{ml/T 2 ) (4.2) 

is a measure of the deviation from conformal symmetry, and this deviation is squared, just as 
discussed in section II A. One could reproduce (4.1) from the derivation of bulk viscosity in 
this paper and of shear viscosity in Ref. [13] by (i) keeping only the |Av 2 | term in the source 
term (3.31), and (ii) making a relaxation-time approximation of the collision operator as a 



A similar estimate was made by Ref. [9] but differs by a factor of 2. The difference is likely due to the 
incorrect identification of the shear viscosity r], by a factor of 2, in Eq. (2.39) of Ref. [9]. 
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rate T(p) that is the same for bulk viscosity and shear viscosity. 17 In Fig. 4, the estimate 
(4.1) is shown by the dashed line for the leading-order result (3.32) for | — v%. It does 
reasonably well at estimating the order of magnitude of our result for bulk viscosity. 

It is interesting that there are certain strongly-coupled but nearly conformal theories 
which find a parametrically different dependence on | — vl than (4.1). In certain theories 
with gravity duals that make them amenable to calculation, 18 Ref. [25] finds £ ~ 77 (| — v£). 
This result is difficult to understand from the picture of viscosity developed in weakly- 
coupled field theories and provides an interesting conceptual puzzle for understanding bulk 
viscosity in strongly-coupled but nearly-conformal theories. 



Acknowledgments 



We would like to thank Larry Yaffe and Sangyong Jeon for useful conversations. We also 
thank Larry Yaffe and Andrei Starinets for goading us into finally doing the calculation. 
This work was supported, in part, by the U.S. Department of Energy under Grant No. DE- 
FG02-97ER41027, by the National Sciences and Engineering Research Council of Canada, 
and by le Fonds Nature et Technologies du Quebec. 



APPENDIX A: ORTHOGONALITY OF S TO THE ENERGY ZERO MODE 



In this appendix, we verify that the source derived in this paper, given by (3.9) and 
(3.22), is orthogonal to the energy- changing zero mode xip) °t E p discussed in section HIE. 
Specifically, we show that (S, E p ) = at the order of our calculation. That is, 

? ,„/^ /0(1±/0) ( P _ ? _„«) £ ^ . (A1) 

This can be checked directly using the QCD-specific formulas of Sec. Ill C, but it is instruc- 
tive to give a more general argument. 

Use dfo = — /o(l ± fo)d(PE p ) and v p = 'VpEp to rewrite the orthogonality condition as 

5>/ wrl E - p - v ^ = v2 ^S w? Epdph - (A2) 

We then need the following two, slightly subtle equilibrium relations, which we will discuss 
below: 

a 

d,e = $>| ^ Ep d p f . (A4) 



Specifically, the ratio of the sources in the two cases then becomes q^/q v — |Ai>g|. The relaxation-time 
approximation is [C\](p) = F{p)x{p)- Using C = 2Qmax\ l=Qi from (3.34) and 77 = j^Q max \ l=2q=qn 
from Ref. [13], one then obtains = I5q^/q^ — 15\Avg\ 2 . 
For other bulk viscosity results in strongly interacting theories, see Ref. [26]. 
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S = Si 



ideal 



The orthogonality relation is then equivalent to the equilibrium relation 

d f3 P-v%e = 0, (A5) 

which is satisfied because v 2 = dP/de = (dpP)/(dpe). 

For the rest of this appendix, we will use the short-hand notation f to stand for 

Deriving general relations for pressure and energy density and their derivatives in a gas 
of quasi-particles is slightly subtle because the effective energies E p of the quasi-particles 
depend on temperature and include the effects of interactions with other quasi-particles. 
The energy density is not simply e = J E p f , for example, because this expression suffers 
the usual Hartree problem of double-counting the interaction energy. [And, if e actually 
were f E p f , we would not get (A4) because there would be an additional term where the 
dp hit the E p .] As discussed in Refs. [21, 22], one simple solution to this problem is to start 
with the entropy density S rather than P or e. Up to higher-order corrections which we 
shall review in a moment, the entropy density of a quasi-particle gas is given by the naive 
ideal gas formula, 

/3(Pidcal + Cideal) = P J {\p • v p + E p ) f = -(3 J \E pP ■ V p / , (A6) 

where the last step follows by integrating the term involving v p = V P E P by parts. Starting 
from this formula for the entropy, we can then use the thermodynamic relation S = drP to 
write dpP = -T 2 S and obtain (A3). 

To get the formula for dpe, it is convenient to first use v p f = =FTV p ln(l ± / ) and 
integrate by parts to rewrite (A6) as 

S = J [±ln(l±/ )+/3£ p /o]. (A7) 

Then use the the thermodynamic relations e = TS — P and dpP = —T 2 S to write 

dpe = dp(TS) -dpP = TdpS. (A8) 

Use of (A7) for S then produces the desired formula (A4). 

It remains only to discuss the approximations that have been used in this analysis. In 
evaluating the entropy, the treatment of the system as an ideal gas of on-shell propagating 
quasi-particles breaks down at order g 3 and above. (See, for instance, the analysis in Ref. 
[22].) But it is adequate to obtain the 0(g 2 ) and the 0{mfy terms in the entropy. For 
massless QCD, that might sound inadequate, because the breaking of conformal invariance 
is an 0(g 4 ) effect. For example, the effective energy of a hard quark is given by 

E 2 p ^p 2 + \ g 2 {T) T 2 =p 2 + I g 2 ^) T 2 + Ift g\p) T 2 ln(T//i) + • • • , (A9) 

and it is the last term which breaks conformal invariance. However, this 0(g 4 ) conformal- 
breaking log is determined by knowledge of the 0(g 2 ) contribution; any 0(g 3 ) or additional 
0(g A ) contributions to thermodynamic quantities will be conformal, up to corrections of 
0(g 5 ), and so will not contribute to the leading-order bulk viscosity. 
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